Frustrated Further-Neighbor Antiferromagnetic and Electron-Hopping Interactions in 
the d = 3 tJ Model: Finite- Temperature Global Phase Diagrams from 

Renormalization-Group Theory 
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The renormalization-group theory of the d — 3 tj model is extended to further-neighbor an- 
tiferromagnetic or electron-hopping interactions, including the ranges of frustration. The global 
phase diagram of each model is calculated for the entire ranges of temperatures, electron densities, 
further/first-neighbor interaction-strength ratios. With the inclusion of further- neighbor interac- 
tions, an extremely rich phase diagram structure is found and is explained by competing and frus- 
trated interactions. In addition to the Ttj phase seen in earlier studies of the nearest-neighbor d = 3 
tj model, the rnt phase seen before in the d — 3 Hubbard model appears both near and away from 
half-filling. 
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I. INTRODUCTION 



The simplest model electron conduction model, includ- 
ing nearest-neighbor hopping on a lattice and on-site 
Coulomb repulsion, is the Hubbard model In the 
limit of very strong on-site Coulomb repulsion, second- 
order perturbation theory on the Hubbard model yields 
the tJ model 0, Q, in which sites doubly occupied by 
electrons do not exist. Studies of the Hubbard model 
4 1 and of the tJ model [sj, including spatial anisotropy 
ij and quenched non-magnetic impurities 0] in good 
agreement with experiments, have shown the effective- 
ness of renormalization-group theory, especially in cal- 
culating phase diagrams at finite temperatures for the 
entire range of electron densities in d = 3. These calcula- 
tions have revealed new phases, dubbed the r phases, 
which occur only in these electronic conduction mod- 
els under doping conditions. The telltale characteristics 
of the T phases are, in contrast to all other phases of 
the systems, a non-zero electron-hopping probability at 
the largest length scales (at the renormalization-group 
thermodynamic-sink fixed points) and the divergence of 
the electron-hopping constant t under repeated rescal- 
ings. Furthermore, the phase diagram topologies, the 
doping ranges, and the contrasting quantitative t and 
antiferromagnetic behaviors under quenched impurities 
7| have been in agreement with experimental findings 
g, 9]. A benchmark for this renormalization-group ap- 
proach has also been established by a detailed and suc- 
cessful comparison, with the exact numerical results of 
the quantum transfer matrix method [lol . [ll| , of the spe- 
cific heat, charge susceptibility, and magnetic suscepti- 
bility in ci = I calculated with our method. [l2| Fur- 
thermore, results with this method have indicated that 
no finite-temperature phase transition occurs in the tJ 



model in d = 1. A phase separation at zero tempera- 
ture has been found in d = 1 in Ref. [l^. Thus, the 
d = 1 tJ appears to have a first-order phase transition 
at zero temperature that disappears as soon as temper- 
ature is raised from zero, as in other d = 1 models such 
as the Isin g an d Blume-Capel models [13, [Hi • A phase 
separation [l6l - [la | occurs in d = 2 for low values of t/J, 
but not for t/J > 0.24. [HI In d = 3, a narrow phase 
separation occurs, as seen in the density - temperature 
phase diagrams below. Two distinct t phases have been 
found in the Hubbard model 4], Tnt and tjj, respectively 
occurring at weak and strong coupling. The calculated 
low-temperature behavior and critical exponent of the 
specific heat [3] have pointed to BCS-like and BEC-like 
behaviors, respectively. Only the Ttj phase was found in 
the tJ model. 

The current work addresses the issue of whether both 
T phases can be found in the tJ model, via the inclusion 
of further- neighbor antiferromagnetic (J2) or further- 
neighbor electron hopping (12) interactions. We find 
that, depending on the temperature and doping level, 
the further-neighbor interactions may compete with the 
further-neighbor effects of the nearest-neighbor interac- 
tions, namely that frustration occurs as a function of 
temperature and doping level. This competition (or rein- 
forcement) between the interactions of successive length 
scales underpins the calculated evolution of the phase 
diagrams. Global phase diagrams are obtained for the 
entire ranges of each type of further-neighbor interac- 
tion. With the inclusion of further-neighbor interactions, 
an extremely rich phase diagram structure is found and 
is explained by competing and frustrated interactions. 
Both THb and Ttj phases are indeed found to occur in 
the tJ model with the inclusion of these further-neighbor 
interactions. Furthermore, distinctive lamellar phase di- 
agram structures of antiferromagnetism interestingly sur- 
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round the r phases in the doped regions. 



II. THE tJ HAMILTONIAN 



rescahng factor 6 = 2. By neglecting the noncommuta- 
tivity of the operators beyond three consecutive lattice 
sites, a trace over all states of even-numbered sites can 
be performed p^[20j. 



On a d-dimensional hypercubic lattice, the tJ model is 
defined by the Hamiltonian 



PH = P 



-t 



E(4 



V 



{ij) 



P, (1) 



where /? = l/fcsT and, with no loss of generality Q, 
t > is used. Here and cja are the creation and 
annihilation operators for an electron with spin a =t or 
J, at lattice site i, obeying anticommutation rules, rii — 
n.i^ + riii are the number operators where riia = cl^Cicr, 

and Si — J2aa' cl^^aa' Cia' is the single-site spin opera- 
tor, with s the vector of Pauli spin matrices. The projec- 
tion operator P = Yiii^ ~ ^ii'^it) projects out all states 
with doubly-occupied sites. The interaction constants 
t, J, V and /i correspond to electron hopping, nearest- 
neighbor antiferromagnetic coupling (J > 0), nearest- 
neighbor electron-electron interaction, and chemical po- 
tential, respectively. From rewriting the tJ Hamiltonian 
as a sum of pair Hamiltonians —/3H{i,j), Eq. ([TJ becomes 



— JSi ■ Sj + Vnijij + ^{ui + n.j) 

m 



P (2) 



where fi = p,/2d. The standard tJ Hamiltonian is a spe- 
cial case of Eq. ([2]) with V/J — 1/4, which stems from 
second-order perturbation theory on the Hubbard model 

Si- 



III. RENORMALIZATION-GROUP 
TRANSFORMATION 

A. d — 1 Recursion Relations 

In d ~ 1, the Hamiltonian of Eq. ^ is 

= ^{-/3i^(^,^-Hl)} . (3) 

i 

A decimation eliminates every other one of the successive 
degrees of freedom arrayed in a linear chain, with the 
partition function being conserved, leading to a length 



Tr p~l^" - Tr pT.d-PHiiA+1)} 

-•-^cvcn*^ — J-i even*- 

even even 
i i 

(4) 

where —/3'H' is the renormalized Hamiltonian. This ap- 
proach, where the two approximate steps labeled with 
~ are in opposite directions, has been successful in the 
detailed solutions of quantum spin [l^-[2^ and electronic 
0-01 systems. The anticommutation rules are correctly 
accounted within the three-site segments, at all succes- 
sive length scales, in the iterations of the renormalization- 
group transformation. 

The algebraic content of the decimation in Eq. Q is 



~l3'H'(i.k) ^ rj,^.^-PH(i,j)~pH(3,k) ^ 



(5) 



where fc are three consecutive sites of the unrenor- 
malized linear chain. The renormalized Hamiltonian is 
given by 



-p'H'{i,k) = P 



- J'Si ■ Sk + V'n.nk + + n^) + G' 



P, 



(6) 



where G' is the additive constant per bond, which is 
always generated in renormalization-group transforma- 
tions, does not affect the flow of the other interaction 
constants, and is necessary in the calculation of expec- 
tation values. The values of the renormalized (primed) 
interaction constants appearing in —/3'H' are given by 
the recursion relations extracted from Eq. ([5]) , which will 
be given here in closed form, while Appendix A details 
the derivation of Eq. ([7]) from Eq. ([5]): 



1 ,„ 74 



t = - in— , J = In — , V = - In — , 

2 72 77 4 7|7| 

m' = M+^1ii(^) , G' = 5'^G + In 71, (7) 



_ 1 i„ 7i7677 
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where 71 = 1 + 2m3/(^) , 



74 = 1 + —U V + -U XJ 



899 I0./3J V n 
3J y At 



76 — 2v X + x/ 



2 3,44, . 1^ M 
77 = — wa; H — V +w/ 

" 3 3 ''V8 2 2 



(8) 



Phase 


Interaction constants at sink 


t 


M 


J 




d (dilute disordered) 





~oo 








D (dense disordered) 





00 








AF 

(antiferromagnetic) 





00 


—00 


—00 

J ^ 4 


TtJ 

(BEC-like superconductor) 


00 
-> 1 


00 


00 

^ -> 2 
n 


— 00 

Z _^ _3 


THb 

(BCS-like superconductor) 


— 00 


CO 


—00 


— CXJ 
,7^4 



TABLE I: Interaction constants at the phase sinks. 



and V ^exp{-J/8 + V/2 + fi/2) , 
a; = exp(3J/8 + y/2 + ^/2) , ■u = exp(/i/2) 

A 



f{A) = cosh V2F+A^ 



Sinh V2t2 + yl2 _ 

(9) 



B. d > 1 Recursion Relations 

The Migdal-Kadanoff renormalization-group proce- 
dure generalizes our transformation to d > 1 through 
a bond- moving step ^1, |^ . Eq. dT]) can be expressed as 
a mapping of interaction constants K — {G, t, J, V, /i} 
onto renormalized interaction constants, K' = R(K). 
The Migdal-KadanofF procedure strengthens by a factor 
of b"^^^ the bonds of linear decimation, to account for 
a bond- moving effect [2^, [l^l- The resulting recursion 
relations for d > 1 are, 



K' = 6''-iR(K), 



(10) 



which explicitly are 
in 

72 



t' - 



ln2:i,J' = fe''-iln2:^,F' = 
77 



■In- 



4 4 

7174 



fo''- V + ^ In , G' = b'^G + b''-' In 71. 



(11) 



This approach has been successfully employed in studies 
of a large variety of quantum mechanical and classical 
{e.g., references in [3]) systems. 



C. Calculation of Phase Diagrams 
and Expectation Values 

The global flows of Eq. ([TU]). controlled by stable and 
unstable fixed points, yield the phase diagrams in tem- 
perature versus chemical potential [1^ : The basin of at- 
traction of each fixed point corresponds to a single ther- 
modynamic phase or to a single type of phase transition. 



according to whether the fixed point is completely stable 
(a phase sink) or unstable. Eigenvalue analysis of the re- 
cursion matrix at an unstable fixed point determines the 
order and critical exponents of the phase transitions at 
the corresponding basin. 

Table I gives the interaction constants t, J, V, fi at the 
tJ model phase sinks. The Ttj and THb phases are the 
only regions where the electron-hopping term t does not 
renormalize to zero at the phase sinks. On the contrary, 
in these phases, t ^ 00 and t — >■ — cxi, respectively. 

To compute temperature versus electron-density (dop- 
ing) phase diagrams, thermodynamic densities are cal- 
culated by summing along entire renormalization-group 
flow trajectories. |23| A density, namely the expectation 
value of an operator in the Hamiltonian, is given by 



1 dlnZ 



(12) 



where is an element of K = {Ka}, Z is the parti- 
tion function, and N is the number of lattice sites. The 
recursion relations for densities are 



fi 



where 



/3a 



(13) 



In terms of the density vector M = {Ma} and the recur- 
sion matrix T = {T^q}, 



/ hd- ^ ^ ^ ^ \ 

/ " dt dJ dv \ 

at' dt' dt' dt' 

dt dJ dV dn 

dJ' dJ' dJ' dJ' 

dt dJ dV dn 

dV dv' av' av' 

at aj av a^i 

n. aji! a^i a^i a^' 

" at aj av au. ' 



(14) 



Eq. ([T3|) simply is 



M = b-'^Ml • T. 



(15) 



At a fixed point, the density vector Mq, = M'^ = M* is 
the left eigenvector, with eigenvalue 6'*, of the fixed-point 
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recursion matrix T* (Table II). For non-fixed-points, it- 
erating Eq. p3|) n times, 

M = fe""''M<^") • T*^"^ • T^""^^ T^^) , (16) 

where, for n large enough, the trajectory arrives as close 
as desired to a completely stable (phase-sink) fixed point 
and M(") ~ M*. The latter density vector M* is the left 
eigenvector of the recursion matrix with eigenvalue . 
When two such density vectors exist, the two branches 
of the phase separation of a first-order phase transition 
are obtained 29, 30] , as illustrated with the phase sepa- 
rations found below. 



Phase sinks 


Expectation values at sink 




{Ui) 






d 














D 





1 





1 


AF 





1 


1 

4 


1 


nj 


2 
3 


2 
3 


1 

4 


1 

3 


THb 


0.664 


0.668 


0.084 


0.336 



TABLE II: Expectation values at the phase sinks. The ex- 
pectation values at a sink epitomize the expectation values 
throughout its corresponding phases, because, as explained 
in Sec. IIIC, the expectation values at the phase sink under- 
pin the calculation of the expectation values throughout the 
corresponding phase which is constituted from the basin of 
attraction of the sink. 



IV. FURTHER-NEIGHBOR INTERACTIONS, 
TEMPERATURE- AND DOPING-DEPENDENT 
FRUSTRATION AND GLOBAL PHASE 
DIAGRAMS 

For the results presented below, we use the theoret- 
ically and experimentally dictated initial conditions of 
V/J = 1/4 and t/J = 2.25. 

The details of the thermodynamic phases found in this 
work, listed in Tables I and II, have been discussedprevi- 
ously within context of the nearest-neighbor tJ 043 3i.nd, 
for the THb phase, Hubbard [3] models. The Tnt phase 
is seen here in the tJ model with the inclusion of the 
further-neighbor antiferromagnetic or electron-hopping 
interaction. Suffice it to recall here that the r phases 
are the only phases in which: (1) the electron- hopping 
strength does not renormalize to zero, but to infinity; 

(2) the electron density does not renormalize to com- 
plete emptiness or complete filling, but to partial empti- 
ness/filling, leaving room for electron/hole conductivity; 

(3) the nearest-neighbor electron occupation probabil- 
ity does not renormalize to zero or unity, again leav- 
ing room for conductivity at the largest length scales; 

(4) the electron-hopping expectation value is non-zero at 
the largest length scales; (5) the experimentally observed 
chemical potential shift as a function of doping occurs 




K 



FIG. 1: Construction of the further-neighbor models. Part of 
a single plane of the three-dimensional model studied here is 
shown. 

^; and (6) a low level (^ 6%) of quenched non-magnetic 
impurities causes total disappearance, in contrast to the 
antiferromagnetic phase (~ 40% for total disappearance) 
[3, again as seen experimentally. The low-temperature 
behavior and critical exponent of the specific heat Q 
have pointed to BCS-like and BEC-like behaviors for the 
THb and Ttj phases, respectively. 

The only approximations in obtaining the results be- 
low are the Suzuki- Takano and Migdal-KadanofF proce- 
dures, explained above in Sees. Ill A and IIIB respec- 
tively. There are no further assumptions in Secs.IVA and 
IVB below. 



A. The t2 Model 

The t2 model includes further-neighbor electron- 
hopping interaction, as shown in Fig. 1. The three- 
site Hamiltonian, between the lattice nodes at the lowest 
length scale, has the form: 

-/3i7(z, j, fc) = - pH{t,j) - fiH{j, k) 

a 

where —/3H{i,j) is given in Eq. so that the first 
equation of Eq. ([7]) gets modified as 

i' = iln^+i2, (18) 
2 72 

only for the first renormalization. Thus, for d — 3, the 
first equation of Eq. (ITT]) gets modified as 

t' = 2 1n— +4t2, (19) 

72 

only for the first renormalization. Thus, the hopping- 
strength t2 contributes to the first renormalization, but 



5 




FIG. 2: Global phase diagram of the further-neighbor t2 model for t/J = 2.25 in temperature vs. chemical potential (first 
column) and, correspondingly, temperature versus electron density (second column). The relation t/J — 2.25 is used for 
all renormalization-group trajectory initial conditions. The t2/t values are given in boxes. The dilute disordered (d), dense 
disordered (D), antiferromagnetic AF (lightly colored), nj (medium colored), and THb (darkly colored) phases are seen. Second- 
order phase transitions are drawn with full curves, first-order transitions with dotted curves. Phase separation occurs between 
the dense (D) and dilute (d) disordered phases, in the unmarked areas within the dotted curves in the electron density vs. 
temperature diagrams. These areas are bounded, on the right and and the left, by the two branches of phase separation densities, 
evaluated by renormalization-group theory as explained in Sec.IIIC. Note that these coexistence regions between dense (D) 
and dilute (d) disordered phases are very narrow. Dashed curves are not phase transitions, but disorder lines between the 
dense and dilute disordered phases. As explained in the text, on each side of the thick full curves (not a phase boundary), 
the further-neighbor electron hopping affects the t phases oppositely. On the dash-dotted curve (also not a phase boundary; 
overlaps, for t2/t = 0, with the thick full curve) electron hopping in the system is frustrated. 



is not regenerated by this first rcnormalization. Note 
that the quantitative memory of the further-neighbor in- 
teraction is kept in all subsequent renormalization-group 
steps, as the flows are modified by the different val- 
ues of the first-renormalized interactions due to the ef- 
fect of the further-neighbor interaction. The subsequent 



global renormalization-group fiows are in the space of 
t, J, V, ^, as is the case in position-space renormalization- 
group treatments j31r.33j using a prefacing transforma- 
tion. Which surfaces in this large (4-dimensional) flow 
space of t, J, V, ^ are accessed is controlled by the original 
further-neighbor interaction. Thus, the further-neighbor 
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FIG. 3: The continuation of the 

interaction t2 shifts the value of t' obtained after the 
first renormahzation-group transformation, as dictated 
by the physical model (Fig.l). Since the value of the 
first-renormalized t' in the absence of t2 already has a 
complicated dependence on the unrenormalized temper- 
ature and electron density, the variety of phase diagrams 
is obtained. Indeed, the effect of the further-neighbor 
interaction is dependent on the electron density, temper- 
ature, and other interactions in the system, due to the 
presence of the first-term in Eq. which is the key 

to the resulting spectacularly different phase diagrams 
as the further-neighbor interaction is varied. (1) If the 
two terms in Eq. (IT5|) are of the same sign, the nearest- 
neighbor and further-neighbor electron hopping terms of 




0.4 0.5 0.6 0.7 0.8 0.9 



<n.> 
I 

global phase diagram in Fig. 2. 

the original system reinforce each other and the r phases 
are enhanced. (2) If the two terms are of opposite signs, 
the nearest-neighbor and further-neighbor electron hop- 
ping terms of the original system compete with each other 
and, with the introduction of further-neighbor electron 
hopping, the r phases are initially suppressed, but en- 
hanced as further- neighbor hopping becomes dominant. 
The two regimes (1) and (2) are separated by the thick 
full lines in the phase diagrams in Figs. 2 and 3. In 
the case (2) of opposite signs, when the two terms cancel 
out each other, the system is frustrated, in which case, 
after the first renormalization, there is no electron hop- 
ping in the system. Since this condition is closed under 
renormalization, both on physical grounds and of course 
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FIG. 4: Global phase diagrams of the further-neighbor J2 model for t/J = 2.25 in temperature vs. chemical potential (first 
column) and, correspondingly, temperature versus electron density (second column). The relation tjj = 2.25 is used for 
all renormalization-group trajectory initial conditions. The J2/J values are given in boxes. The dilute disordered (d), dense 
disordered (D), antiferromagnetic AF (lightly colored), rtj (medium colored), and rnh (darkly colored) phases are seen. Second- 
order phase transitions are drawn with full curves, first-order transitions with dotted curves. Phase separation occurs between 
the dense (D) and dilute (d) disordered phases, in the unmarked areas within the dotted curves in the electron density vs. 
temperature diagrams. These areas are bounded, on the right and and the left, by the two branches of phase separation densities, 
evaluated by renormalization-group theory as explained in Sec.IIIC. Note that these coexistence regions between dense (D) 
and dilute (d) disordered phases are very narrow. Dashed curves are not phase transitions, but disorder lines between the 
dense and dilute disordered phases. As explained in the text, on each side of the thick full curves (not a phase boundary), 
the further-neighbor interaction affects the antiferromagnetic phase oppositely. On the dash-dotted curve (also not a phase 
boundary; overlaps, for J2/J = 0, with the thick full curve), antiferromagnetism in the system is frustrated. 



in our recursion relations (Eq. (O), no t phase can oc- 
cur in such a system. The dash-dotted curves in Figs. 2 
and 3 indeed show such systems. These competition, re- 
inforcement, and frustration effects are temperature and 



doping dependent. These, and all other physical effects, 
do not depend on the sign of nearest-neighbor t of the 
original unrenormalized system, due to the symmetry of 
hypercubic lattices [Sj] and as seen in Eq. 
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2.5 0.5 




FIG. 5: The continuation of the global phase diagrams in Fig. 4. 



Figs. 2 and 3 give the global phase diagram of the 
t2 model, as a function of temperature, electron den- 
sity, chemical potential, and t2/t. The values of the 
hopping-strength ratios t2/t for the consecutive panels in 
these figures are chosen so that they sequentially produce 
the qualitatively different phase-diagram cross-sections, 
thereby revealing the evolution in the global phase dia- 
gram. Second-order phase transitions are drawn with full 
curves, first-order transitions with dotted curves. Phase 
separation occurs between the dense (D) and dilute (d) 
disordered phases, in the unmarked areas within the dot- 
ted curves in the electron density vs. temperature dia- 
grams. These areas are bounded, on the right and and 



the left, by the two branches of phase separation den- 
sities, evaluated by renormalization-group theory as ex- 
plained in Sec. inc. Note that these coexistence regions 
between dense (D) and dilute (d) disordered phases are 
very narrow. 

The cross-section ^2 = is the phase diagram ob- 
tained in previous work [5]. This diagram contains the 
Ttj phase between 33 — 37% hole doping away from half- 
filling and at temperatures 1/t < 0.12. The thick full 
curve here gives the systems devoid of electron hopping 
due to the combined effects of temperature and doping 
on a nearest-neighbor-only interaction system. The first 
term of Eq. is positive on the high density/chemical 
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potential, low temperature side of the thick full curve 
and negative on the low chemical potential/density, high 
temperature side of the thick full curve. Thus, the in- 
clusion of ^2 > will create competition and frustra- 
tion (respectively reducing and eliminating the r phases) 
on the low chemical potential/density, high temperature 
side of the curve discussed here, reinforcement (enhanc- 
ing the T phases) on the high chemical potential/density, 
low temperature side of the same curve. The opposite 
occurs at i2 < 0. The thick full (no- hopping) curve of 
t2 = is included, again as thick and full, in the i2 7^ 
phase diagrams and the effects discussed here are seen 
in the evolution, in both directions, of the global phase 
diagram. 

Pursuing the negative values of ^2, we see at t2/t — 
—0.0625 that the Ttj phase, being below the thick full 
curve, is indeed reduced and bisected into two discon- 
nected regions by the frustration (dash-dotted) curve. At 
the more negative value of t2/t = -0.125, only the higher 
doping region of the tjj phase remains and is enhanced 
as explained after Eq. (ITS)) , extending through a wider 
range to 45 — 55% hole doping. The antifcrromagnetic 
and disordered phases take part in a complex lamellar 
structure, in a narrow band between 35 — 45% hole dop- 
ing at low temperatures. At the even more negative val- 
ues of t2lt = —0.25 and —0.5, the Ttj phase appears in 
a wide range of hole doping, between 35 — 55%. Besides 
the complex lamellar structure of antifcrromagnetic and 
disordered phases, we also see that the Tnb phase par- 
ticipates in the lamellar phase structure and, separately, 
appears adjacently to the antiferromagnetic phase near 
half-filling. Particularly near half-filling, the THb phase 
which evolves adjacently to the antiferromagnetic phase 
reaches to the higher temperatures of 1/t ~ 0.5. This 
topology is identical to that obtained for the Hubbard 
model [4]. 

For the positive values of t2/t, the r phases are en- 
hanced as explained after Eq. ^5)) and the topology 
quickly evolves to that encountered in the Hubbard 
model. The Ttj is not bisected by the frustration (dash- 
dotted) curve and appears between 33 — 37% hole dop- 
ing as a continuation of the structure at t2 = 0. The 
Tub phase occurs again in two distinct regions and the 
one which lies nearer to half-filling again extends to high 
temperatures. 



B. The J2 Model 

The J2 model includes further-neighbor antiferromag- 
netic interaction, as shown in Fig. 1. The three-site 
Hamiltonian, between the lattice nodes at the lowest 
length scale, has the form: 

-PH{^,].k) ^- pH{i,3)- pH{3,k) 

-J2^S,.S., (20) 



where —f3H{i,j) is given in Eq. ([2]), so that the second 
equation of Eq. ([7]) gets modified as 

J' = ln^ + J2, (21) 

77 

only for the first renormalization. Thus, for d = 3, the 
second equation of Eq. ([TT]) gets modified as 

J' = 41n^+4J2, (22) 
77 

only for the first renormalization. Again, the interaction 
J2 contributes to the first renormalization, but is not 
regenerated by this first renormalization. Reinforcement 
or competition occurs when J2 is, respectively, of same 
or opposite sign as the first term in Eq. (22). These two 
regimes are again separated by the thick full lines in the 
phase diagrams of Figs. 3 and 4, while again frustration 
occurs on the dash-dotted lines. In the reinforcement 
regime, we expect a large extent of the antiferromagnetic 
phase. The Tnt phase is also expected to grow in the 
reinforced region, for it is found along the temperature 
extent of the antiferromagnetic phase. 

Figs. 4 and 5 show the global phase diagram of the 
J2 model, as a function of temperature, electron den- 
sity, chemical potential, and J2/J. Again, the values 
of the coupling-strength ratios J2/J for the consecutive 
panels in these figures are chosen so that they sequen- 
tially produce the qualitatively different phase-diagram 
cross-sections, thereby revealing the evolution in the 
global phase diagram. Again, the phase separation re- 
gions of the first-order phase transitions are very nar- 
row. For negative values of J2/J, the antiferromagnetic 
phase is enhanced, both near half-filling by the mech- 
anism explained after Eq. ([^ cind, separately and to 
a lesser extent, displacing the Ttj phase. The latter 
behavior is similar to that seen under the introduction 
of quenched impurities, both experimentally [S^-HB] and 
from renormalization-group theory [7|. The Tub phase 
improves near the large antiferromagnetic region near 
half-filhng. At J2/J — —2, the th& phase is found in 
a wide range of hole doping, namely between 15 — 30%. 
Another interesting result is that the r^j phase is de- 
pressed in temperature but remains stable in the interval 
of 33 — 37% hole doping. 

For positive values of J2/J, the antiferromagnetic 
phase is reduced in the region near half-filling and en- 
hanced in the region near the Ttj phase, for reasons 
explained after Eq. ([^ . The r^b phase grows adja- 
cently to the enhanced antiferromagnetic region, being 
located above the Ttj phase, causing a complex structure 
at higher hole dopings and low temperatures. 

C. Conclusion 

We have shown that the tJ model with further- 
neighbor antiferromagnetic (J2) or further- neighbor elec- 
tron hopping (12) interactions exhibits extremely rich 



global phase diagrams. The phase separation regions of 
the first-order phase transions are very narrow. Further- 
more, these calculated phase diagrams are understood in 
terms of the competition and frustration of nearest- and 
further-neighbor interactions. We find that the two types 
of r phases, previously seen in the Hubbard model, oc- 
cur in the tJ model with the inclusion of further-neighbor 
interactions. 
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TABLE III: The two-site basis states, with the corresponding 
particle number (n), parity (p), total spin (s), and total spin 
z-component (ms) quantum numbers. The states \(j>3}, l^s), 
and \4)s) are obtained by spin reversal from 102), 104), and 
1 07), respectively. 



Appendix A: Derivation of the Decimation Relations 

The derivation of Eq. ([T]) , first done in Ref. ^] , is given 
in this Appendix. In Eq. ([5]) the operators —[3'H'{i,k) 
and —f3H{i,j) — (3H{j,k) act on two-site and three-site 
states, respectively, where at each site an electron may 
be either with spin a or 4, or may not exist (0 state). 
In terms of matrix elements, 

J2{y'^w,Vk\e-^''^'^^^-^''^J^''^\u,w,v,) , (Al) 

Wj 

where Ui,Wj,Vk,Ui,Vf, are single-site state variables, so 
that the left-hand side reflects a 9 x 9 and the right-hand 
side a 27 X 27 matrix. Basis states that are simultaneous 
eigenstates of total particle number (n), parity (p), to- 
tal spin magnitude (s), and total spin z-component (mg) 
block-diagonalize Eq. (jAip and thereby make it man- 
ageable. These sets of 9 two-site and 27 three-site eigen- 
states, denoted by {It/ip)} and {|V'g)} respectively, are 
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TABLE IV: The three-site basis states, with the correspond- 
ing particle number (n), parity (p), total spin (s), and total 
spin 2-component {iris) quantum numbers. The states |'04-5), 
10)7), 1015-16), 1019), |0'2i), 10)23), 1026-27) are obtained by 
spin reversal from \tp2-3), {tpe), |0ii-i2), \tpn), IV'20), 1^022), 
|'024-25), respectively. 
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TABLE V: Block-diagonal matrix of the renormalized two-site 
Hamiltonian —l3'H'{i,k). The Hamiltonian being invariant 
under spin-reversal, the spin-flipped matrix elements are not 
shown. 



given in Tables III and IV. Eq. (|A1[) is thus rewritten as 

(0p|e-'''^'(-'=)|0p) = 

J2 E('^pl"^^'^)("»^J-«fel^9)(^g|e"''''^''^'^"''''^^'''^lV',-)- 

{^q\u,WjVk){UiVk\(l)p) ■ (A2) 

There are five independent elements for 
(0p|e-'^'^'(*''=)|(/)p) in Eq.dJll (thereby leading to five 
renormalized interaction constants {t', J', T^', /x', G"}), 
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which we label 7p, 

7p^(0p|e-^'^'(^''=)|</.p) forp= 1,2,4,6,7. (A3) 

The diagonal matrix {(j)p\—P'H'{i, k)\(l)p) is given in Table 
V. The exponential of this matrix yields the five renor- 
malized interaction constants in terms of "fp , as given in 
Eq. ([71) ■ Furthermore, according to Eq. (jA2p . each 7^ is 
a linear combination of some {'ipq\e~'^-'^^^'^^~'^-'^^^'''^tpq), 

71 = (V'l||V'l) + (V'2||l^2) + (V'4||V'4), 

72 = (V'3||^3) + ^('08||V'8) + (l/'l2||V'l2) + ^(V'l3||V'l3) , 

74 = (V'6||'06) + ^('09||'09) + (l/'l7||V'l7) + ^('018||'018) , 

76 = (V'io||'0io)+2(V'22||V'22> , 

2 4 

77 = (V'll||'01l) + 2(^2o||V'2o)+g('024||V'24) , 



TABLE VI: Diagonal matrix blocks of the unrenormalized 
three-site Hamiltonian ~l3H{i,j) — l3H{j,k). The Hamilto- 
nian being invariant under spin-reversal, the spin-flipped ma- 
trix elements are not shown. 



where (V'qll^/'g) = {^Jq\e-^"^'-^^-^"^^-''^ijq) . In order to 
calculate (^/'qle"'^^'''-'^"'^^'-''''^ the matrix blocks in 
Table VI are numerically exponentiated. 
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